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We give a general definition for the tunneling time in the Landau-Zener model. This definition 
allows us to compute numerically the Landau-Zener tunneling time at any sweeping rate without 
ambiguity. We have also obtained analytical results in both the adiabatic limit and the sudden limit. 
Whenever applicable, our results are compared to previous results and they are in good agreement. 
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I. INTRODUCTION 

Tunneling is one of many fundamental quantum pro- 
cesses that have no classical counterparts. It exists ubiq- 
uitously in quantum systems and is a key to under- 
standing many quantum phenomena [l|, |3, H, 0i S S 
0, H, H, lis El, Discussions on tunneling can be 

found in all textbooks on quantum mechanics. How- 
ever, these discussions are mainly focused on the prob- 
ability of tunneling from one quantum state to another 
or from one side of a potential barrier to the other. In 
contrast, there are only a few extensive and in-depth 
discussions in literature on another aspect of tunneling, 
the time of tunnelin g, t hat is, how long a tunneling pro- 
cess lasts [H, [il [H, III [Ullil. This disparity is partly 
caused by the difficulty to define properly tunneling times 
in many situations. 

The difficulty of having a proper definition for tunnel- 
ing time has its root in the wave or probabilistic nature 
of quantum mechanics, and is best demonstrated in the 
example of a wave packet tunneling through a potential 
barrier. In this case, one would intuitively define the 
tunneling time as the time spent by the peak (or cen- 
troid) of the wave packet under the barrier. However, 
as pointed out by Landauer and Martin [Toj. if this defi- 
nition was used, a packet could leave the barrier before 
entering it. To overcome this difficulty, many different 
definitions have been suggested, and no clear consensus 
has been reached so far[2(j|. 

The focus of this study is on the tunneling time in the 
Landau-Zener (LZ) model[2l|, [l^. This system is much 
simpler than the wave packet and barrier system. Nev- 
ertheless, a proper definition of tunneling time in this 
model is still missing in spite of the studies in the past 
by many authors. Mullen et al. [1^ discussed the LZ 
tunneling time in the diabatic basis for the two limit- 
ing cases, the adabatic limit and the sudden limit. They 
found that for large A (the adiabatic limit), the tunnel- 
ing time scales with A/a and for small A (the sudden 
limit), the tunneling time is about ^Jhja. A is the min- 
imal energy gap between the two eigenstates in the LZ 
model and a is the sweeping rate. However, Mullen et al. 
did not giv e a general definition for the tunneling time. 
Vitanov|24| has given a much more thorough study on 
the LZ tunneling time. He obtained analytical results 



for the LZ tunneling time in both diabatic and adiabatic 
bases. Furthermore, Vitanov tried to give a general def- 
inition for the tunneling time. However, his definition 
fails in certain cases, in particular, the adiabatic limit. 

In this paper, we give a general definition for the tun- 
neling time in the LZ model. We show both analytically 
and numerically that this definition yields reasonable re- 
sults at any sweeping rate in both adiabatic basis and 
diabatic basis. With this general definition, we are able 
to reproduce the previous results obtained by Mullen et 
aZ.;23| and Vitanov 24 1. Furthermore, we are able to com- 
pute analytically the tunneling time in the adiabatic basis 
at the adiabatic limit, which is given by 2\/2i/3 _ lA/a. 
This, to our best knowledge, has not been obtained be- 
fore. 



Besides it theoretical significance, our work has also 
potential applications. In the Monte-Carlo simulation of 
quantum tunneling in molecular magnets, the authors in 
Ref. [25| have used an empirical formula for the LZ tun- 
neling time. This formula, given by ^ A^/a^ -I- 2h/a and 
interpolating the two limiting results in ReLfSsI, is not 
well founded. With this newly-proposed definition, one 
no longer needs this empirical formula to do the Monte- 
Carlo simulation. 



We note here that it is important to study the tun- 
neling time in both the diabatic basis and the adiabatic 
basis. When the LZ model is applied to the case of the 
tunneling of spin under a sweeping magnetic field [1^, the 
diabatic basis is a better choice. When the LZ model is 
applied to the tunneling between Bloch bands under a 
constant force, it is better to use the adiabatic basis [26j. 

Our paper is organized as follows. In Sec|TTl we intro- 
duce our definition of the tunneling time in the LZ model 
and we analyze the effectiveness of our definition. In Sec. 
mil we present our results of the tunneling times, which 
include the analytical results at the adiabatic limit and 
the sudden limit and the numerical results for the gen- 
eral case. Our results are given both in the diabatic basis 
and the adiabatic basis. In the last section, we discuss 
our results and conclude. 
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FIG. 1: (a)The solid line is the time evolution of the probabil- 
ity function Pd (t) and the dashed line is its step-like function 
fit. (b) The time derivative of the two functions in (a) and 
T] = 0.2565. 




FIG. 2: Time evolution of tunneling probability Pa{t) and its 
time derivative dPa{t)/dt at the adiabatic limit, (a) The time 
evolution of the tunneling probability Pa{t)\ (b) the function 

of dPa{t)/dt, T) = 0.10. 
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II. DEFINITION OF THE TUNNELING TIME 

The LZ model is a two-level system and is described 
by [2ll,[2| 



dt \b{t) 



a{t) 
b{t) 



where 



(1) 



(2) 



with 7 = ckf changing with time linearly, a is usually 
called sweeping rate. The tunneling behavior in the LZ 
model can be described in two different bases, diabatic 
basis and adiabatic basisjl^l. In the diabatic basis, we 
use Pd{t) = \b{t)\'^ to describe the LZ tunneling dynam- 
ics. In the adiabatic basis, we use 



Pa{t)^\a*{t)a{j{t)) + b*mjm' 



(3) 



where (0(7), 6(7)) is the instantaneous eigenstate of 
H{^). When our discussion is independent of the ba- 
sis, we will remove the subscript and simply use P for 
both Pd and Pa- 

The time evolution of the probability function P{t) can 
be found by numerically solving Eq.(IT]). A typical result 
of P{t) is shown in Fig[Tl where we see a sharp transition 
occurs around t = and is followed by decaying oscilla- 
tions. This observation suggests an intuitive (or natural) 
definition for the LZ tunneling time. One may first fit 
the P curve with a smooth step-like function (dashed 
line in Fig[T]Ja), and then define the half- width of its 
time derivative (dashed line in Figlljb)) as the tunneling 
time. However, like its counterpart in the wave packet 
and barrier system, this intuitive definition of tunneling 
time fails because of its two shortcomings. First, there 
are numerous methods to find the fitting step-like func- 
tion in Fig[lja); there is no obvious criterion by which 



one method is better than the other. Secondly, at the adi- 
abatic limit, the P curve looks drastically different from 
the typical case in FiglIJa). As shown in Figl^l the P is a 
single-peaked function; one really has to be far-stretched 
to fit it with a step-like function. These two drawbacks 
show that this intuitive definition of the tunneling time 
based curve-fitting is not a good choice. One has to find 
an alternative. 

In Ref . |24i| , Vitanov introduced a general definition for 
the LZ tunneling time. His definition is 



P(oo) 



(4) 



where the time derivative value P (0) at t = is used to 
represent the rate of the transition around t = 0. For the 
typical time evolution of the probability function P(t) 
shown in Fig[Tl this definition works well. However, Vi- 
tanov's definition fails at the adiabatic limit like the in- 
tuitive definition as we shall see later. In fact, Vitanov's 
definition does not work in the adiabatic basis in general. 
A typical P curve in the adiabatic basis is shown in FiglSl 
where P'(0) clearly over-represents the transition rate. 

We overcome these difficulties and find a general def- 
inition of the tunneling time for the LZ model. In the 
definition, we first find a i' < such that 



P(^ ) — 2 "^^^^2; ; 



(5) 



where P„ 



is the maximum value of P when t < 0. 



Usually, Pmax = P{t = 0). The above condition will 
be called half-width condition from now on for ease of 
reference. We then introduce two more variables 



/■° d 

Si=J ■^^Pit)dt = P(0) 



^2 



dt 



P{t)dt = P{oo) - P(0) , 



(6) 



(7) 



3 



0.15 



"ra 0.1 




0.05 



FIG. 3: Time evolution of the tunneling probability in the 
adiabatic basis, rj = 1.425. 



which are the left {t < 0) area and right {t > 0) area of 
the dP/dt curve, respectively. With these defined vari- 
ables, we define the tunneling time as 



T=\t'\{l 



(8) 



Three quick remarks, i) The three variables in the def- 
inition can be computed without any ambiguity, ii) For 
the typical case shown in Fig[Tl we have Si ~ S2 and, 
therefore, r = 2|i'|, which is in agreement of the "intu- 
itive" definition that we discussed before, in) The abso- 
lute value is used because S2 can be negative in certain 
cases, for example, the case in Fig 15] The reason of the 
appearance of negative S2 is that the transition around 
t — over-shoots the overall transition P{oo) — P{~(X)) 
and the system needs to spend some time to "wind back" . 

In the following, we shall apply our definition and com- 
pute the tunneling times in the LZ model. Both analyti- 
cal and numerical approaches will be used. The analyti- 
cal approach is used for two limiting cases, the adiabatic 
limit and the sudden limit. For the LZ model, we can 
introduce a "quickness" parameter 



2ha 



(9) 



The adiabatic limit is 77 ^ 1 while 77 ^ 1 corresponds to 
the sudden limit. For a general case, we have to resort 
to the numerical method. We first solve numerically the 
equation of motion Eq.([T]), then compute the tunneling 
probability function P{t), and finally find the tunneling 
time with our definition in Eq.(l8]). In our computation, 
we use A = 1.2x 10~^fcb, where kb is Boltzmaniiconstant, 
which is a typical value in molecular magnet Jl 23. 

The tunneling time will be computed in both the adi- 
abatic basis and the diabatic basis. For clarity, we shall 



use Td for the tunneling time in the diabatic basis and 
for the tunneling time in the adiabatic basis. 
III. TUNNELING TIMES IN THE DIABATIC 
BASIS 

We first consider the diabatic basis and follow it with 
the discussion on the adiabatic basis in the next section. 



A. Analytical results 

At the adiabatic limit [r] ^ 1), according to 
VitanovfH 



1 



at 



2 2Va2i2 + A2 ' 



(10) 



the variable t' can be obtained from the half-width con- 
dition 



Pd{t)\t=t' =\{Pd{t <G)Uax ^\ 



The result is 



t' 



\/3 A 
3 a 



(11) 



(12) 



Since we have Si — S2 for the tunneling curve Eq. pil)) . 
the tunneling time with our definition of Eq. ^ is 



a_ 2V3A 



3 a 



(13) 



which agrees well with the result of Mullen et al. [23| . 

At the sudden limit rj 1, it is beneficial to take a 
transformation 



a<2 

a{t) = a(<)exp(-i — ) , 



(14) 
(15) 



for the LZ model. As a result, the diagonal terms in the 
Hamiltonian are transformed away and we can expand 
a{t) and b{t) in powers of ry (effectively. A) For the 
initial condition a(— 00) = 1 and b{—oo) = 0, we obtain 

°° 1 2 

ait) = [1 + E(-l)' (^«2fe(2/)] exp(-z^) , (16) 

bit) = [E(-l)'^'^(^^2.-i(.)] cxp(z^) , (17) 



fe=i 



where 
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an{y) / exp{ixl/2) / exp(-ixl/2) , 



K{y) = / exp(-ia;i/2) / ex]i{ixl/2) . 



exp[(— l)"^"'^ia;^/2](ia;i(ix2 . . . dxn-idxn 



exp[(— l)"ia;„/2](ia;ic?X2 . . . dxn-idxr, 



(18) 
(19) 



with y = t/ {h/aY^^ . At the sudden limit, it is sufficient 
to keep Eq. ([17]) to the lowest order of 1/ ry. Consequently, 
we obtain 



_ 1 

~ 2r] 



1 



exp(— i — )(ix 



(20) 



exp(i — ) + / exp(— I — )da; 

where Ciyl-Jir) and S{y/^/Tr) are the Fresnel 
integrals [231 • One can prove that the maximum 
value of Pd for t < is at t = 0, where Pd ~ ^/{Ari). 
Thus with the half-width condition 



Pdit') = \PdiO)^^^, 



(21) 



we find numerically that t' ~ — 0.6241y ^. According to 
Refs.im,!!!, we have 



Pd{oo) = 1 - exp( ) . 

Therefore, at the sudden limit ( 77 3> 1), we have 



Pd(oo) - Prf(O) 



Based on our definition in Eq. ([| 



(22) 



(23) 



A\t'\ = 2.4964a/- 
a 



the tunneling time is 



,rf^.|.|-....u.,,-, (24) 
V a 

which agrees well with the result of Mullen et al. [l^l ■ 



B. Numerical results 

Our numerical results of the tunneling time in the 
diabatic basis is plotted in the log-log scale in Fig(3| In 
the figure, we see that the results for the two limiting 
cases are connected by a smooth kink. We have also 
compared these results with the empirical relation ~ 
\J 1^ jd^ + 2h/a that was used in Ref . [1^ ; the agreement 
is quite good. 

We have amplified the results at the adiabatic and the 
sudden limits and plotted them in Fig[5| and Fig[6l re- 
spectively. In these two figures, we have also compared 
them to the analytical results and the agreement is ex- 
cellent. 




FIG. 4: Tunneling times r^j in the diabatic basis. The solid 
line is our numerical results and the dashed line is the empir- 
ical formula — A^/q^ -|- 2h/a used in Ref. [2^. a is in 

the unit of =r. 
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FIG. 5: Tunneling time at the adiabatic limit in the di- 
abatic basis. The circles are the theoretical results given in 
Ea. (|13|l and the solid line is the numerical results, a is in the 



unit of 



IV. TUNNELING TIMES IN THE ADIABATIC 
BASIS 



As mentioned already, when one applies the LZ model 
to describe the tunneling between Bloch bands, it is more 
convenient to use the adiabatic basis. It turns out that 
the results in the adiabatic basis are quite different from 
the ones in the diabatic basis. 
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FIG. 6: Tunneling time at the sudden limit in the diabatic 
basis. The circles are the theoretical results given in Ea. (l24p 
and the solid line is the numerical results, a is in the unit of 

2fi ■ 



A. Analytical results 

As in the case of the diabatic basis, at the adiabatic 
Umit(?7 <C 1), the tunneling probability function Pa in 
the adiabatic basis has been found by Vitanov 



AiaH"^ + A2)3 ■ 
With the half-width condition, 



Pa{t)\t=t' ^l{Pa{t<0)) 



we find that 



t' = -V2V3_l- 
a 



max 



A 



(25) 



(26) 



(27) 



Thus, based on our definition in Eq.®, the tunneling 
time is 



t: = 2v/2i/3 _ i_ , 



(28) 



which is very similar to the tunneling time Eg. psp in 
the diabatic basis. In contrast, according to Vitanov's 
definition Eq.Q, the tunneling time is^24| 



V2A , ttA\ 



/ah 



iah' 



(29) 



At the adiabatic hmit (a 0), the tunneling time 
tends to be zero. This result contradicts with the physi- 
cal reality, the tunneling time should be very long at the 
adiabatic limit. Thus, Vitanov's definition does not work 
in this case. Alternatively, our result in Eq. can be 
viewed as the first successful attempt to find the tunnel- 
ing time at the adiabatic limit in the adiabatic basis. 

We next consider the sudden limit (77 1). In terms of 
y = t/ y^h/a, the instantaneous eigenstates of the Hamil- 
tonian ^ are 




1/2N 

y Ml/2 



(30) 



According to Eq. (fTB]) and Eq. pT|) . we can obtain the 
tunneling probability up to the first order of l/?7 



Pa{y) 

With Eq. (fT9| ■ we arrive at 



ci(y) c,{y)) ( "^P^-^^^l 

\hi{y) eyi-p{i\) ^ 



(31) 



2 v/'^lw' + 2) 



cos— / sm— da; + —(- — - 

2 2 27]4 2 jy2 



y 



(32) 



It is quite obvious that the last two terms of the above 
equation is much smaller than the first two terms. Fi- 
nally, we obtain 

(33) 



^' 2 2 Va2t2 + A2 



which is surprisingly identical to the result at the adia- 
batic limit in the diabatic basis (see Eq. fTU)) '). As a result, 
we can similarly obtain the tunneling time 

3 a 

which agrees very well with Vitanov's result ^2^ . 



B. Numerical results 

In the adiabatic basis, the numerical results of the tun- 
neling time are shown in Fig|71 We see that the results in 
the two limiting cases are also connected by a kink. How- 
ever, this kink is not as smooth as the kink in the diabatic 
basis; the first derivative of the tunneling time with re- 
spect to a is not continuous. The numerical results for 
the two limits, the adiabatic limit and the sudden limit, 
are plotted and compared to the theoretical results in 
Figl5]and Fig[51 respectively. Again, we find an excellent 
agreement. 
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FIG. 8: Tunneling time Ta at the adiabatic limit in the adi- 
abatic basis. The circles are the theoretical results given by 
Ea. (l28|l : the solid line is the numerical results, a is in the 

A 2 

unit of =r-. 
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FIG. 9: Tunneling time Ta at the sudden limit in the adiabatic 
basis. The circles are the theoretical results given by Ea. (|34[) : 
the solid line is the numerical results, a is in the unit of 



V. DISCUSSION AND CONCLUSION 

In the above, we have obtained analytical results for 
the tunneling times in the LZ model at two different lim- 
its and in two different bases. They are, respectively, f^, 
t% and t^. We have found that only is proportional 
to yjhja while the rest of the three tunneling times all 
scale as A/a. It is not hard to understand why t^, the 
tunneling time at the sudden limit and in the diabatic 
basis, does not scale as A/a. The effect of A is to couple 
the two bare states, (1,0) and (0, 1). which serve as the 
base vectors in the diabatic basis. At the sudden limit, 
the system changes very fast and its wave function re- 
mains almost unchanged. As a result, the system does 
not feel the effect of A. It is also not hard to understand 
that the two tunneling times at the adiabatic limit, iJJ 
and i^, scale as A/a. At the adiabatic limit, the effect 
of A is fully felt by system and gets reflected in the tun- 
neling time. 

The most puzzling is f*, the tunneling time at the 
sudden limit in the adiabatic basis. Unlike the other 
tunneling time at the sudden limit, it is proportional 
to A/a. Moreover, its corresponding probability func- 
tion Pa{t) described by Eg. ([55)1 is surprisingly identical 
to the probability function Pd{t) in Eg. lfTO)) . which is at 
the adiabatic limit in the diabatic basis. To understand 
this, we have to look into the details of the evolution. 
At the adiabatic limit, the system follows its instanta- 
neous eigenstate as demanded by the quantum adiabatic 
theorem pSj. Pd{t) in Eq. ([TO|) is obtained by projecting 
this instantaneous eigenstate to the bare state (0, 1). At 
the sudden limit, the wave function of the system changes 
little and remains in the bare state (1,0). However, in 
the adiabatic basis, this wave function needs to be pro- 
jected to the instantaneous eigenstate to obtain Pa{t) de- 
scribed by Eg. ([33)1 . As we know, projecting a bare state 
to an instantaneous eigenstate is the identical to project- 
ing the same instantaneous eigenstate to the same bare 
state. This explains why the probability function Pa{t) 
in Eo. ([33)) is the same as Pd{t) in Eg. pO]) . Consequently, 
this also explains why scales as A /a. 

In sum, we have presented a general definition of the 
tunneling time for the Landau-Zener model. We have 
shown that this definition works for any sweeping rate 
and can be used for the numerical computation of the 
tunneling time without any ambiguity. In particular, 
we have obtained analytical results for the two limiting 
cases, the adiabatic limit and the sudden limit. We have 
not only reproduced known results but also found the 
tunneling time at the adiabatic limit in the adiabatic ba- 
sis, which has not been found before to our knowledge. 
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